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Abstract 

We consider here a layered superconductor subject to an externally applied moderately-strong 
electromagnetic field. We predict hysteretic jumps in the dependence of the surface reactance of 
the superconductor on the amplitude Hq of the incident electromagnetic wave. This very unusual 
nonlinear phenomenon can be observed in thin superconducting slabs at not very strong ac ampli- 
tudes, if the frequency of the irradiating field is close to the Josephson plasma frequency. Using the 
set of coupled sine-Gordon equations, we derive the expression for the phase shift x °f the reflected 
wave and obtain the conditions for the appearance of hysteresis in the x(-£^o)-dependence. 

PACS numbers: 74.78.Fk 74.50. +r 
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I. INTRODUCTION 



There has been a recent surge of studies of electromagnetic waves (EMWs) propagating in 
artificially fabricated media (see, e.g., Refs. , including metals with modulated proper- 
ties^, arrays of coupled waveguides 5 , left-hand materials^^, and layered superconductors^. 
The excitation of these waves can produce a large variety of resonance anomalies^ in the 
reflectivity, transmissivity, and absorptivity, offering new types of optical nano-devices. 

The recent increase of these type of studies is related to nonlinear surface and waveguide 
EM modes (see, e.g., Refs. 1 3 7 ) . In this broad context, a layered superconductor is a 
medium favoring the propagation of nonlinear— and surface^ waves in the (important 
for applications^^ 5 -) terahertz (THz) and sub-THz frequency ranges. Both the existence 
of surface waves^ and the nonlinear effects^^ 1 ^ occur due to the gap structure (see, e.g., 



Ref. Il6l ) of the spectrum of Josephson plasma waves, which was experimentally observed 
via Josephson plasma resonance 17 . The nonlinearity of Josephson plasma waves with fre- 
quency u close to the Josephson plasma frequency ujj becomes important even at small field 
amplitudes oc |1 — uj 2 /uj 2 ] 1 ^ 2 . In close analogy to nonlinear optics^, the nonlinear JPWs ex- 
hibit numerous remarkable features^ 1 ^ 1 ^, including the slowing down of light, self- focusing 
effects, and the pumping of weaker waves by stronger ones. However, the nonlinearity of 
EMWs in layered superconductors is quite different from optical nonlinearities. This leads 
one to expect very unusual phenomena in the EMW propagation in this nonlinear media. 

In this paper, we predict and analyze theoretically one of such unexpected nonlinear 
effects in a thin slab of a layered superconductor subject to an externally applied electro- 
magnetic wave. We show that, under specific conditions, the amplitude dependence of the 
phase x °f the reflected wave becomes many-valued. This should result in hysteretic jumps 
of x when sweeping the amplitude of the incident wave, a phenomenon which is very unusual 
for conductors and superconductors. 



II. PROBLEM STATEMENT AND EQUATIONS FOR THE ELECTROMAG- 
NETIC FIELD 

Consider a slab of a layered superconductor of thickness d (see Fig. [Q. The crystal- 
lographic ab-plane coincides with the xy-plane and the c-axis is along the z-axis. The 



2 



interlayer distance D is much smaller than the thickness d of the slab. 
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FIG. 1: (Color online) Geometry of the problem. The slab of a layered superconductor is irradiated 
with a p-polarized electromagnetic wave with the magnetic field symmetrical with respect to the 
middle of the sample. 



Let the sample be irradiated by two p-polarized (Transverse Magnetic) plane monochro- 
matic electromagnetic waves with the magnetic fields symmetric with respect to the mid- 
dle of the sample, the plane z — 0. Therefore, the magnetic H = {0, H, 0} and electric 
E = {E x , 0, E z } fields satisfy the symmetry conditions, 

H(x, z, t) = H(x, —z, i), 
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E x (x, z, t) = -E x (x, -z, t), E z (x, z, t) = E z (x, -z, t). (1) 

Due to this symmetry, we will only consider the field distribution in the upper half-space 
z > 0. 

The electromagnetic field in the vacuum, z > d/2, is the sum of the incident and specu- 
larly reflected waves. The Maxwell equations give for them, 

H v (x,z,t) =if cos7_ + # r cos(7+ + x), 

E^(x,z,t) = - -^# sin7_ - F r sin(7 + + x)), (2) 
7_ = k x x—k z z — cut, 7+ = k x x + k z z — k z d — tut, 

with k x = ksm9, k z = kcos9, k = uj/c. Here u is the wave frequency, c is the speed of 
light, and 9 is the angle of incidence. The value of x is the phase shift of the reflected wave 
at the boundary z — d/2 of the slab. As is known, x defines the surface reactance X of a 
sample: when neglecting the dissipation, 

X = — tan (—} cos 9. 
c \2J 

Inside a layered superconductor, the electromagnetic field is determined by the inter- 
layer gauge-invariant phase difference ip of the order parameter. The spatial distribution of 
<p(x, z, t) obeys the set of coupled sine-Gordon equations (see, e.g., Ref. Il9l]20l ). We consider 
the nonlinear JPWs with \ip\ 1, when the Josephson current J c sunp can be approximated 
by J c ((p — (p 3 /G). We also assume that the gauge-invariant phase difference experiences small 
changes on the scale D in the z-direction, and thus we can use the continuum approach. In 
the continuum limit, the coupled sine-Gordon equation has the form, 

Here X a b and A c = c/uje 1 ^ 2 are the London penetration depths across and along layers, 
respectively, and uj = (87ieDJ c /heY^ 2 is the Josephson plasma frequency. The latter is 
determined by the maximum Josephson current J c , the interlayer dielectric constant e, and 
the interlayer spacing D. The spatial variations in the z-direction of the fields inside the 
very thin superconducting layers are neglected. Here we also omit the dissipation terms 
related to the quasiparticle conductivity. They are controlled by the sample temperature T 
and can be reduced to negligibly small values. Moreover, in Eq. ([3]), we neglect the term 
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with the capacitive coupling, for waves with sufficiently high k x = ksinO ~ uj/c ^> (3/\ c . 
Here (3 = R 2 D e/ sD <d 1 is the prefactor of the capacitive coupling 2 ^, Ru is the Debye length 
for a charge in a superconductor, and s is the thickness of the superconducting layers. 

The magnetic and electric fields in a layered superconductor are related to the gauge- 
invariant phase difference as 



dH s Hp f 1 dV _^ 

F s_ K>&H a _ $ U) 



c dzdt 1 u 2ttDA_ 
where $o — nch/e is the flux quantum and e is the elementary charge. 



As was shown in Ref. |T3|, the nonlinearity in Eq. ([3]) can play a crucial role in the wave 
propagation for frequencies close to uj, i.e., for |1 — Q 2 \ = |1 — u 2 /oj 2 \ <C 1. Indeed, if 
ip ~ |1 — Q 2 \ <C 1, the cubic term <p 3 in Eq. ([3]) is of the same order as the linear term 
tufd 2 <p/dt 2 + (p. 

We consider the frequency range below the Josephson plasma frequency, Q < 1, and seek 
a solution of Eq. ([3]) of the form, 

(p(x, z, t) =a(z)(l - Vl 2 ) 1/2 sin( 7 o + a), 
7o =k x x — k z d/2 — ut, 

keeping only the first harmonics in (k x x — uot) . 

Substituting ip in Eq. (jSJ) into Eq. one obtains 



(5) 



(I _ Q2\ 

H s (x, C, t) = -Ho- '-K0 cos( 7o + «), 

K 

E s x (x, (,t)=Ho il ~ n) P h'(C) sin( 7 o + a). 

K 



(6) 



Here we introduce the dimensionless variables, 



*(C)=a(C)-nR C = f, (7) 



A a 5 ft Ac/Ca; 



and parameters, 

P — " «■ — /-Q\ 

A cV /T (l-fi 2 )V2' 

and the prime denotes d/d(. 

Equations ([3]) and (jSJ) yield the ordinary second-order differential equation for a(z) 

d 2 



1 K \e 



a 3 



— + K 2 a = 0. (9) 



Further, we also assume k ^> 1, which is valid for not very small incident angles 9. In this 
case, integrating Eq. (Q with the symmetry condition, 

a'(0) = 0, 

we obtain 

iM 2 ^^) 2 - 1 ' (10) 

where clq = a(0). The solution of Eq. ffTQ]) can be written in the implicit form, 




The phase diagram, i.e., the set of a\a) curves for different values of the constant do, is 
shown in Fig. [2j Solid circles mark the sample boundaries, while open circles indicate the 
middle of the slab. Arrows show the direction of motion along the phase trajectories when 
increasing the coordinate (. 

Thus the electromagnetic fields in the vacuum and in the superconducting slab are de- 
termined by Eqs. (j2J) and (jSJ), respectively. The latter equations contain the function a(Q 
given by Eq. (flTl) . 



III. PHASE SHIFT OF THE REFLECTED WAVE 

Now we find the relationship between the phase shift x °f the reflected wave and the 
amplitude H of the incident wave. With this purpose, we join the tangential components 
of the electric and magnetic fields in the vacuum and in the superconducting slab, at the 
interface z = d/2. Thus, separating terms with sm(k x x — ut) and cos(k x x — out), we derive 
four equations for H r , if (d/2), x, and a: 

—h + h r cos x —Ph'(8) sin a, 
—h Q — h r cos x =h{5) cos a, 

(12) 

—hr sin x =Ph'(5) cos a, 
—hr sin x —h(6) sin a. 

Here 

Hn K H„ K 




FIG. 2: (Color online) The phase diagram a' (a) (for a > 0). Recall that a is the amplitude of the 
gauge-invariant phase ip (see Eq. ([5])), while a' is its derivative with respect to the dimensionless 
coordinate Q = KzjX^. The paths along the arrows on the phase trajectories between solid circles 
correspond to the change of the coordinate z inside the sample, from z = —d/2 to z = d/2. Open 
circles correspond to the middle of the slab (z = 0). 



5 = nd/2\ a i ) is the value of ( at the interface z = d/2. Note that, without any loss of 
generality, we can assume the value h(S) of the total magnetic field at the sample surface to 
be positive, h(S) > 0. A negative h(5) correspond to replacing \ ~ ¥ X + tt in Eqs. ( IT2l) . 



7 



Excluding h r and a from Eqs. (fT2|) . we find 



o = \\lh?{5) + {Ph>(5)) 2 



(13) 



X = 2 arctan 




(14) 



Equations (TT31 . ffT41) . together with Eqs. ([7j) and (TTTi) . give, in an implicit form, the required 
dependence of the phase shift x on the amplitude of the incident wave Hq. 

For further analysis of the amplitude dependence of the phase shift Xi it is very important 
to take into account the non-single- valued relation between the values h(S) and a(S). Indeed, 
Eq. ([7]) and Fig. [3] show that there exist three values of a that correspond to the same value 
of h, if < h < (32/27) l ' 2 . Taking into account Eq. ( TT3l) . we conclude that three different 
values of the parameter h'(5)/h(5) in Eq. (TBI) correspond to the same value of h^. This 
results in the appearance of three branches of the dependence x{ho)- 

Below, for simplicity, we restrict ourselves to the case of small sample thicknesses, when 



and one can easily derive the asymptotic equations for all three branches x(^o)- 
A. Low-amplitude branch x(^o) 

First, we discuss the branch of the dependence x(^-o) that corresponds to the portion 
"1-2" of the h(a) curve in Fig. 3. This monotonically increasing branch is shown by the 
blue curve (with number 1) in the main panel of Fig. 4, and also in inset (a) of Fig. 4. This 
blue branch is defined within the interval (0, ho cr ) of the ho change. The ending point ho c ? 
of the branch is described by the expression, 



<J< 1. 



(15) 



In this case, \a(() — ao\ 1, and Eq. (jTTj) is significantly simplified, 




(16) 




(17) 



At this point 




(18) 
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FIG. 3: (Color online) The normalized magnetic field h versus the normalized amplitude a of the 
gauge-invariant phase (p. This h{a) dependence is described by Eq. ([7]). Arrows indicate the motion 
along the h(a) curve when changing the normalized amplitude ho of the incident wave. 



The low-amplitude asymptotics of x(^o) dependence is, 

X-2Ps(l + ^Y /i «l. (19) 
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FIG. 4: (Color online) The numerically-obtained dependence of the phase shift x °f the reflected 
wave on the dimensionless amplitude tiQ of the incident electromagnetic wave. The values of the 
parameters here are: d/X^ = 0.05, A c /A a b = 200, k = 10, and e = 16. The insets show magnified 
portions of the x(^o) plots. The arrows on the curves indicate how the phase shift x changes 
when periodically varying the amplitude Hq. The first jump in x(ho) a t ho = /locr is shown by 
the downwards vertical arrow in the inset (b). This produces a jump ~ —9P5. The second 
(reverse) jump at ho = /i™ m is shown in the main frame and in the inst (c). Namely, the upward 
vertical arrow (from x(M I before jump = to x(Mlaft er jump ~ °) shows the J um P in x(.h )- 
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B. High-amplitude branch %(^o) 



The second and third branches of the x(ho) dependence correspond to the portions "2-3" 
and "4-5", respectively, of the h(a) curve in Fig. 3. The second and third branches of x(ho) 
are shown by the green and red curves (indicated by numbers 2 and 3, respectively) in Fig. 4. 
The second branch exists in the interval 

h™ n <h < h 0cT , h™ n ~ V2P5 . (20) 

As is seen in Fig. 4, the first and second branches of x(ho) almost meet at the point h 0cr . 
The third branch is defined for 

h™ n <h <oo. (21) 

Near the minimum value of h , i.e., for h ~ h™ in , the second and third branches have 
the asymptotics, 

X(ho) =±(*- ^-4P^ . (22) 

Here, the signs "+" and "— " correspond to the second and third branches, respectively. 
The third branch tends to zero for h — > oo following the expression, 

X (h ) = -P5^y /3 , /i »l. (23) 

Now we present the parametrically-defined formula for all three branches, which is valid 
for ho not very close to h 0cr : 

ho(ao) = ^=^1, (24a) 
16P 

X(ao) = o ~ 2 S- (24b) 

8 - a 

Figure 5 demonstrates the very good agreement of this formula with numerical results. 
C. Hysteresis jumps in the dependence x(ho) 

From the analysis shown above, we can now describe the behavior of the phase shift X of 
the reflected wave when periodically changing the amplitude ho of the incident wave. 

When increasing ho from zero, the phase shift x{ho) increases monotonically following 
the first branch in Fig. 4. At the point h = h 0cr , the first branch comes to the end, and a 
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FIG. 5: (Color online) The phase shift x °f the reflected wave versus the dimensionless amplitude 
ho of the incident wave, for d/X a b = 0.05, A c /A a f, = 200, k = 10, and e = 16. The solid curves are 
the plots of the three branches of the function given by Eqs. (|24|) . while the dots were numerically 
obtained from Eqs. (|7l). ([TT]) . (fT3|) . and ([Til). 



jump to the third branch should occur when further increasing ho. The phase shift x(^-o) 
jumps from the value 

x(Ziocr-O) ~ 3P5 

to the final value 

x (h 0ci + 0) ~ -GPS. 
12 



Thus, the jump here is 

AxL = , 0cr - -9PS. (25) 

When decreasing ho, the phase shift x(^-o) decreases following the third branch, crosses 
the point h = h 0cT , and only at the point ho = /i™ m it performs a reverse jump to the first 
branch. Therefore, now 

A X L 0= ^n « tt. (26) 

Thus, the hysteretic jumps A% of the phase x(^o) of the reflected wave could be observed 
when periodically changing the amplitude ho of the incident wave. 



IV. CONCLUSION 



In this paper, we have predicted and theoretically analyzed unusual phenomenon for 
conducting media. We have shown that, due to the specific nonlinearity of layered super- 
conductors, hysteretic jumps of the surface reactance could be observed when periodically 
changing the amplitude H of the incident wave. A remarkable feature of the predicted 
phenomenon is the relatively small values of the necessary ac amplitudes Ho- The hysteretic 
jumps can be observed even when ^< 1. According to our analysis, the critical amplitude 

Hocr = h cT ■ H - *C H ~ 20 Oe. 

K 

There are two small parameter here, (1 — Q 2 ) and 1/k. As was shown in Ref. y, the heating 
effect is negligible for such ac amplitudes. 

The phenomenon discussed here is another exciting example of the numerous unusual 
effects related to the very specific nonlinearity of layered superconductors. 
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